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1 INTRODUCTION 



The random walk model for neutral polymer is perhaps the simplest mathematical rep- 
resentation for long flexible chains |T|]. The tremendous progress made in the theoretical 
understanding of conformational and dynamical properties of flexible polymer systems be- 
comes possible because systematic calculations using the random walk model in conjunction 
with the inclusion of excluded volume interactions can be carried out at least in principle 0]. 
The resulting model, referred to as the Edwards model, is the minimal representation of real 
polymers that adequately describes the global properties of several polymeric systems. The 
random walk model views the flexible polymer chains as a Brownian curve. In the discrete 
representation, a flexible chain can be modeled as one for which angles between successive 
chain segments are not correlated. Since the orientations of chain segments are independent, 
the segment vectors have the Markovian property so that the mean squared end-to-end dis- 
tance is proportional to the number N of segments of size a in the chain. In the continuum 
limit this chain becomes a Brownian curve. 

In contrast there are a number of issues concerning the behavior of semiflexible chains 
that are not satisfactorily solved. Many polymeric molecules have internal stiffness and hence 
can not modeled as freely jointed chains This is especially the case in several biopoljmaers 
such as actin, DNA, and microtubles 0]. A measure of the stiffness of a polymer is in terms of 
the so called persistence length, Ip, which gives an estimate of the length scale over which the 
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tangent vectors along the contour of the chain backbone are correlated. The typical values 
of Ip of biopolymers lie in the range of a several nm to few mm. This range is several orders 
of magnitude larger than the persistence length of flexible chains. If the contour length, L, 
of such molecules is of the same order of magnitude as Ip then it is imperative to include 
bending rigidity to describe the conformations of the chain. In these situations the chain 
behaves as wormlike and the appropriate model for this was introduced sometime ago |^. 
It has been known for a long time that the wormlike model provides a good starting point 
in the theoretical description of chains with internal stiffness. Inspired by several recent 
experiments which have probed a number of properties of stiff biological molecules, there 
has been a renewed interest in understanding their shapes. The purpose of this chapter is to 
provide a simple way of calculating a number of interesting properties of semiflexible chains 
using a simple meanfield variational approach. 

The effect of excluded volume profoundly changes the properties of flexible chains. Al- 
though historically the importance of excluded volume was recognized sometime ago, the 
introduction of the Edwards to represent the effect of this short range interaction made pos- 
sible systematic calculation of various static and dynamic properties using field theoretical 
methods In this chapter, in which we focus on the properties of semiflexible chains, 
the excluded volume effects will be ignored. This is physically reasonable for relatively stiff 
chains for which deviations from rod-like conformations are negligible. 

A simple way to account for the stiffness of a semiflexible chain is to constrain the angles 
between two successive segments 6 to be fixed. The value of 6 depends on the local stiffness 
of the chain. This prescription leads to the freely rotating chain model. If we describe 
the configurations of a polymer chain by the set of position vectors {r„} = (fq, ...,rjv) or 
alternatively by the set of segment vectors {Ar„} = (ri — tq, fat — tn-i), then the spatial 
correlation, (Ar„ ■ Ar„_i), in the freely rotating chain has the assigned value a^cos 6. In 
the continuous limit {a ^ 0,6 ^ 0,N ^ oo, Na = L) the freely rotating chain becomes 
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the so called wormlike chain 0]. In this case, the ratio 2a/ 9"^ defines the persistence length 
Ip, which is the typical length scale over which the chain changes its direction appreciably. 
Other conformational properties of such a model are well known in the literature [3,6-20]. 

The spatial correlations (Ar„ ■ Ar,„), which characterize the properties of a semifiexible 
chain, decay exponentially as exp[— a|?2 — m|//p]. Thus the conformational properties of 
a semifiexible chain beyond the length scale Ip reduce to those of fiexible chains, i.e., one 
can view the stiff chain as being made up of several rigid segments of length Ip that are 
freely joined. However because of the intrinsic skeletal stiffness of many synthetic polymers 
as well as biopolymers one needs to develop a model that explicitly builds effects due to 
chain bending. The chain stiffness turns out to be a relevant parameter in the description of 



isotropic- nematic transition condition in liquid-crystalline polymers |jT5[ . Even for an isolated 
chain, the chain stiffness should be taken into account to describe the local properties of stiff 
polymer chains. This is especially important in polyelectrolytes. The scaling behavior of the 



electrostatic persistence length 1^ is known to depend on the rigidity of the chain [|T^, |T7|. 
Many biological molecules and short chains of otherwise flexible chains also belong to the 
class for which the chain stiffness plays an important role. 

A number of theoretical models have been introduced in the literature to account for 
chain stiffness. The earliest model for stiff chains is the wormlike chain (also known as 
Kratky-Porod model) in which the angles between successive chains are constrained p|. 
Although physically reasonable, this model has not yielded analytically tractable results for 
equilibrium and dynamical properties. Harris and Hearst introduced a "simplifled model" 
of stiff chains in which the tangent vector u(s) = dr/ds was allowed to fluctuate as opposed 
to having the constraint u^(s) = 1 for all s [0]. It has been noted that the resulting model 
does not satisfy the spatial homogeneity of stiff chains. More recently a model that does 
not suffer from this restriction was proposed by Lagowski, Noolandi, and Nickel |jl8| using 
a functional integral formalism. These authors showed that the resulting model yielded the 
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mean squared end-to-end distance in agreement with Kratky-Porod. The spatial correlations 
decay exponentially with a slightly shorter value of the persistence length. 

In this chapter we show that a model for stiff chains proposed by Lagowski, Noolandi and 
Nickel (LNN) |jl8| results from a stationary phase evaluation of certain functional integrals 
that occur in an appropriate field theory for stiff chains. Our approach is systematic and 
can be applied to many diverse problems involving semiflexible chains. We should note that 
Winkler et al. [|19| have also obtained a model for stiff chains using the maximum entropy 
principle. These authors did not notice that their model in the continuum limit is identical 
to that of LNN. Furthermore their method appears more cumbersome than the standard 
functional integral approach presented here. 

The rest of the chapter is organized as follows. In section 2, we present the basic strategy 
behind the meanfield variational approach. The resulting model, as mentioned earlier, leads 
to the LNN representation of wormlike chains. The theoretical ideas are used to calculate 
the distribution of end-to-end distance in semiflexible chains in section 3. The application of 
the theory to the problem of semiflexible chains under tension with forces or the stretching 
of DNA is developed in section 4. The possible limitations of the theory are illustrated in 
section 5 by studying the behavior of semiflexible chains in a stretching nematic field. The 
chapter is concluded in section 6 with a few additional remarks. 



2 MEANFIELD MODEL 
2.1 flexible chains 

The basic methodology can be illustrated using the simpler example of a flexible chain. 
This is the limiting case of a stiff chain as the rigidity vanishes. The probability function for 
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the flexible chain conformations without excluded volume interactions can be written as 



N 



n=l 



where ip = 5(|Ar| — a)/47ra^ denotes the random distribution of a segment vector of length 
a. Eq. (1) accounts for the chain connectivity. We can rewrite the probability weight in Eq. 
(1) by introducing auxiliary fields A„ as 

^ A, 



rioo N r ^ A 

vl>{r„}oc / n^Anexp - ^ ^((Ar, 

•^-«°o„=l L n=l ^ 



(2) 



We now show that a stationary phase evaluation of the free energy of the chain described 
by above weight leads to the probability weight for the Brownian chain. This approximation 
amounts to relaxing the locally enforced constraint of (Ar„)^ = to a global one, ((Ar„)^) = 

,and the validity of the approximation can be justified a posteriori as a — > 0. The free 
energy F of a non-interacting flexible chain can be written as 



Y\ dA„exp(— JF{A„}) 

n=l 

where the free energy functional ^{A„} is defined by 

-^{A„} = -ln<^ / Yl dr^exp -a~^ ^ A„r^ \ -a^^K. 

[J j^^i I J J — 1 

N 

(|lnA„ — Xnoj + const. 



(3) 



n=l 



(4) 



n=l 



In the above equation, the order of the r„ and A„ integrations is interchanged. If we denote 
the trajectory A„ along which the integrand in Eq. has its maximum value by X^, then the 
free energy can be expanded around this stationary phase trajectory, A^'. In the following the 
subscript 'cl' will be omitted. In the meanfield theory for which the constraint is imposed 
only on an average, we retain only the leading term in this expansion. By setting the 
partial derivative of the free energy functional J^{A„} with respect to A„ to zero, we get the 
stationary phase condition: 



The independence of A„ on n reflects the symmetry of the problem of an ideal flexible chain. 
Since the delta function can be also represented as 



the stationary phase evaluation becomes very accurate in the continuum limit, a — 0. 
Thus long flexible chains, i.e., ^ 1 can be well described by the following weight in the 
continuum limit 



where \E'[r(s)] is written in the functional integral notation and is the Wiener measure. 

By treating the random fields A(s) at the mean field level, the microscopic constraints 
conjugate to the fields A(s), which ensure that the chain segments are connected but oth- 
erwise randomly distributed, are relaxed to the global ones. This results in the expected 
probability weight given in Eq. (|^ for a long flexible chain and is the Wiener measure ob- 
tained in the path integral description of a diffusion equation. 

2.2 Linear stiff chains 

The approach described above can be extended to semiflexible chains. In this calculation 
we assume that the stretching of two connected chain segments is not important so that the 
coupling between this degree of freedom and the bending degree of freedom can be ignored 
[^. In this case, the weight in Eq. (|l|) needs to be modified so that it yields non- vanishing 
correlations (Ar„ ■ Ar„_i) = = 2a^/lp. This can be achieved if we multiply the weight 
in Eq. (|l|) by the Boltzmann weight exp(/pa~'^ J2n=i Ar„+i ■ Ar„) corresponding to the local 
interactions between adjacent segments. This term favors parallel alignment of adjacent 
segments over bent configurations. In the A„ representation of the probability weight, this 
can be rewritten as exp[— |/pa~^(Ar„_|_i — Ar„)^] with a redefinition of A^. Then the weight 




(6) 




(7) 
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associated with a particular configuration of a semiflexible cliain becomes 

rioo N r ^ A / 

v[/{r4 oc / n dA„ exp - — ((Ar„)' - «') - ^Ta E (^r.+i - Ar, 

■>-^^n=l L n=l n=l 

In tlie continuum limit, this can be written as functional integral 



^[u(s)] cx exp 



'2" 



L n 5(u^(^)-i) (9) 



0<s<L 

where u(s) = dr{s)/ds is a unit tangent vector. The properties associated with the weight 
\E'[u(s)] are well known in the literature [2-3,10]. The random variable u(s) describes the 
rotational Brownian motion on a unit sphere, = 1. If we let P(us, u^/; s', s) be the 
probability that u(s') = u^' when u(s) = u^, then this function obeys a diffusion equation 
on the unit sphere. The solution of the diffusion equation can be expanded in terms of 
spherical harmonics. This enables us to compute the following correlation 

(u(.')-u(s))=exp(-|s'-s|/g. (10) 

This correlation along with the Markovian property of u leads to the mean squared end-to- 
end distance given by 

{R^) = /^rdsds'(u(s') •u(s)) 
Jo Jo 

= 2/pL-2/p2(l -e-^/'p). (11) 

Even though the results given in Eq. (p!oD and Eq. ([TTD are exact, the use of Eq. (P) to 
describe non-ideal semiflexible chains turns out to be quite formidable. The major difficulty 
arises because of the constraint u^(s) = 1. One encounters similar difficulty in other physical 
systems described by the non-linear cr model [^] for which the magnitude of a spin S is held 



fixed, = const. Thus it is of practical interest to obtain a tractable model for such 
constrained systems. We will extend the stationary phase approach adopted for the flexible 
chain to obtain a tractable meanfield model for a semiflexible chain. 

In our stationary phase approach, the field A„ is treated as a parameter to be determined. 
The dependence of A„ on n depends on the problem under consideration. The free energy 



functional for an ideal semiflexible chain can be written as 

. N \ E ^ 

^{K] = -In / Jl dr„ exp -^-^ + « 

where in a matrix form is given by 

Ea 



(12) 



(13) 



kBT 

with C = {ri, rAr}"^. The 3A^ x 3A^ matrix Q is defined by 



Qnm — ^nSn m — 2^^^ ~^ ^" mil)- (14) 



Then the free energy JF is given by 
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JF{A„} = - ln(det Q) - aJ2K + const. (15) 

2 n=l 

The stationary phase evaluation of amounts to minimizing the free energy with respect 



to A„, i.e 



d , „ 3 91n(detQ) ^ ,^ ^, 

-J^{X^} = ^ ;^ = a, 1 < n < A^. (16) 



It can be easily shown that the minimization condition in Eq. (|1^) amounts to requiring 
(u^) = 1 in the continuous limit. This follows because Eq. ( [TB| ) can be rewritten as ^|-JF = 
((Ar„)^/a^) — 1. This is a simultaneous equation for the unknown parameters A„ for which 
we can not find an analytical solution. An examination of the structure of the matrix Q, 
however, leads to the following properties of A„ which satisfy the above equation; Ai = Aat 7^ 
A2 =,...,= Aat-i. For our purposes it suffices if A„ can be chosen so that (u^(s)) = 1 and 
other conformational properties are reproduced. If all A„ are equal to each other, as is the 
case for the flexible chain, then the chain described by the probability weight in Eq. (^ 
shows inhomogeneity, i.e., the chain fluctuates more strongly at both ends than elsewhere. 
Having recognized the translational asymmetry in the problem of a semiflexible chain, it is 
convenient to rewrite An as follows; Ai = Aa? = A + 5/a, A„ = A (2 < n < — 1). With 
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these simplifications, the weight for the semifiexible chain at the level of a stationary phase 
approximation, becomes 



^mf[u(s)] cx exp -A dsu^(s) - '^^\'ds^ ~ '^('^o + "i) 



;i7) 



The functional in Eq. (0) is identical in form to that proposed by LNN. The explicit 
expression for det Q and thus the stationary-point conditions for A and 5 can be obtained by 
setting a recursion relation in A^. Alternatively, we can exploit an analogy between the path 
integral in Eq. ([T7| ) and the harmonic oscillator in quantum mechanics ||2^ . The propagator 



of a harmonic oscillator of a mass /p and a frequency Vt = Y'2A//p, denoted by Z{uq, ul, L), 
is given by 

[uj^ + Uq) cosh QL — 2uq ■ ul' 



, /27r sinh ^]L^3/2 
Z(uo, u^; L) = — j exp 



Qlp/2 ■ sinh QL 



We can thus rewrite the free energy as 

T[X, 5] = -In J duodui e-^("o+"' ) Z(uo, u^; L) - (LA + 25) + const 



fin 



6 sinh L, 



2A XL , ^ 



2A^ 









I XL 



^ ^""V^" ^ Inj^sinh 



2A^ 



(LA + 25) + const 



(19) 



where we have used 



r r dx dj/e-p'("'+^')±''"^ 



(20) 



A little algebra leads to the following stationary phase conditions for A and 5 which require 
1-^ = = i^: 



/ A^ ^ 3 
2 4 



(21) 



Note here that the values of A and 5 don't depend on the contour length L of a chain. Since 
the stationary phase condition is imposed on A, only one of /p or A is independent. 
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To understand the features implied by the weight in Eq. (|T3), let us compute the corre- 
lation (u(s) ■ u(s')). Using the Markovian property of u, the correlation can be computed: 

(u(s) ■ u(s')) = Af-^ J V[u] u(s) ■ u(s') ^mf[u] 

= AT-i J duodu,du,,dUie-^("o+"i)Z(ui, u,,; L - s') 

xZ(u^/,u^; s' - s)u(s) ■ u(s')^(us,Uo; s) 

= ^ \n\f du,u,,e-^("'+"'')+""-"='Z(u,,, u,; s' - s) 
da a=o U J 

= exp(-|s'-s|//o) (22) 

where A/" is the normalization constant given by A/" = / T'[u] \E'mf[u] and Iq = |/p. A direct 
consequence of the above correlation is (u^) = 1 and thus the constraint = 1 is enforced 
only on an average in the meanfield model of semiflexible chains with the weight given in 
Eq. ([TtI). The above correlation can be compared with one in Eq. ( p!o|) obtained with the 
exact weight. A comparison of Eq. (0) and Eq. (p2D shows that the persistence length for 
the approximate model for stiff chains (cf. Eq. (|T^) is smaller by a factor of |. The plausible 
reason for this is the following: In the original model the constraint condition is u^(s) = 1 
for all values of s. The model obtained by enforcing the global condition (u^) = 1 allows for 
unrestricted (restricted only on an average) fluctuations in u thus allowing for configurations 
that would be prohibited by the restricted condition u^(s) = 1. Thus we would expect Iq to 
be less than /p. 

Winkler et al. have obtained exactly the same result (see their Eq. (4.18)) using the 
maximum entropy principle. Since they also only enforce the constraint on an average the 
resulting theory, as described here, should be viewed as meanfield hke. If Iq is understood 
as a new definition of the persistence length, then the stationary phase weight in Eq. (|T^ 
predicts the same conformational behavior as the exact one in Eq. (^. If dimension d is 
large enough, we expect the effect of fluctuations around the mean value to become mini- 
mal. It is known that the stationary phase approach described here becomes exact as c? — > oo 
[18]. It's worth noting that for a chain described by this weight with 6 = 0, the above corre- 

10 



lation holds only for <^ s, s' ^ L. This is because of the excess end fluctuations in this case. 



2.3 Closed stiff chains 



For practical purposes it's more convenient to use a translationally symmetric model 
for a semifiexible chain, i.e., one described by \E'mf[u(s)] in Eq. (|17D with 6 = 0. For ring 
polymers for which the periodic condition u(0) = u(L) and the closure relation Jq u(s)ds = 
are imposed, we expect all the A(s) to be equal or A(s) should be independent of s. The free 
energy in this case can be written as 



^[A] = -In / P[u(s)]5(/o^ds u(s))^mf[u(s) 



5=0 



-LA. 



(23) 



To compute this integral, it is convenient to introduce the Fourier transform 

.271372' 



Completeness of this expansion is reflected in the closure condition 



(24) 



— J ds exp 



.27rs 

> 

' L 



(n — n] 



6nn.' • 



(25) 



In terms of new variables u„, the weight can be written as 



^mf{u„} = exp<j - J2 



2 L 



+ XL 



Ur, ■ U 



(26) 



where the contribution from n = is excluded to incorporate the closure condition /q^ u(s)ds = 0. 
Consequently free energy becomes 



^[A] = -In / n du„|det|£^| ^mf{u„} - LA. 



Since the transformation in Eq. (EJ) linear, the Jacobi determinant 



det^ 



(27) 



is independent 



of u„. Furthermore it does not depend on A. Its value is thus unimportant and will give rise 
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to an additive constant to the free energy: 



^[A] 



—In j du„\E'MF{u„} — LX + const. 



XL' 



= 3 In ^sinhL^ 



+ const. 



A 



- I In(LA) - LA + const. 



(28) 



where const, refers to A-independent terms. In the last step of the above equation, we have 
used 

In sinh X = In x + ^Inh + ). (29) 

n=l ^ ^ ""^ 

Now the stationary phase condition reads 



'^ = |(cothL. 



A 



1 21 



\2l. 




(30) 



Here we have L-dependent condition for A. This is because only paths which satisfy the 
cychc conditions, u(L) = u(0) and /Q^u(s)ds, contribute to the free energy J^. The cychc 
conditions are also incorporated in the correlation of u(s), i.e., 

.271 



(u(s') • u(s)) = {un- u„/)exp 



i — {s'n' + sn) 
Ij 



IT.- 



exp 



_|_ lp(2nnf 
2 L 



(31) 



where (...) is defined by 

Upon using the following identity 



E 



cos nx 



I Un^O dUnUn ' Un'^MpfUn} 
/ rin^O dUn*MF{Un} 



TT cosh a{Tr — x) 1 



(32) 



+ 2a sinh ctTr 2a^ 
we can rewrite the correlation of the unit tangent vectors in a closed form 



(33) 



(u(5') • n{s)) 



cosh[r)(I - 2|,s' - ,s|)/2] - 2/r)I • siiih(r)I/2) 
cosh(OL/2) - 2/VtL ■ sinh(0L/2) 



(34) 
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Note here that (u(s + L) ■ u(s)) = (u^(s)) = 1. In the hmit of L — oo, however, the cydic 
conditions are irrelevant. That is the stationary phase condition and the correlation given 
above reduce to those of open chains. This can be checked by taking the limit L oo in 
above equations for the ring polymers. An alternative method (without derivation) for ring 
segments has been proposed earlier [^. This involves modifying the original Harris- Hearst 
model for open chains. Huber et al. have computed quasi-elastic scattering for a modified 



version of this model 24 



2.4 Linear periodic chains 

The ring polymer described by Eq. ( ^Bf ) was introduced to circumvent difficulties as- 
sociated with inhomogeneity in the linear chain with a uniform stationary phase value A. 
According to the correlation function in Eq. ([31|), however, we can assume a uniform A even 
for the chains with periodic boundary conditions only without violating the homogeneity of 
the chain. In this case we should include n = contribution in the Fourier modes; the free 
energy is thus given by 



^[A] = 3 InJ^sinhL^ 
This leads to the following stationary phase condition 



) - LA + const. (35) 



\/¥-t"thI^^. (36) 
The periodic boundary condition is also incorporated in the correlation function which reads 

(u(.).u(/))^ -''l'^-JI;'-f»/^l (37) 

cosh[\lL / 2) 

This correlation shows that the chain with periodic boundary condition is homogeneous; 
(u^(s)) = 1 for all s, < s < L. As L —* oo, the periodic boundary condition is irrelevant 
as in the ring polymers. 
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An alternative and useful presentation of this calculation can be made writing the cur- 
vature term in a symmetric way. Due to the periodic boundary conditions, i.e., u(0) = u(L) 
and f u(0) = ^u(L), the free energy functional for the semiflexible chain can be written as 

J^[A] = -ln J Vu{s)exp £ dsds'u{s)Q{s,s')u{s) - J ds\{s) (38) 



where the operator Q is defined by 

Q{s,s') = 



Sis'-s). 



(39) 



Integration with respect to u(s) leads to 



J^[X] = I tr In dsA(s). 



(40) 



The stationary phase condition can be obtained by requiring dT /d\{s) = 0: 

2n ,/ 1 \ 



(41) 



2 + '^('^')/ s,s 

If a uniform stationary phase value A is assumed, the above equation can be Fourier trans- 
formed. To this end, it is convenient to define a complete set of eigenstates {\s)} with s 
a curvilinear space label, and a complete set Fourier-conjugate to this, such that 

^|n) = i^^\n) . Thus \s) and \n) are related with each other through {n\s) = -^exp[i^^]. 
With the aid of these, the right hand side of the above equation can be written as 



" 2 



I E ( 

n=— oo 
oo 

1 E 



1 



\n) {n\s) 



\ r I ip I 
n=— oo •^-i-' ' 2 ' L 



(42) 



(27rn)2 • 
2 L 

Combined with the identity in Eq. (|33D , the stationary phase condition relation can be 
rewritten as 

(43) 



= I cothL. 



A_ 
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This condition becomes the same one as in Eq. (pT]). The model described by Eq. (35) with 
L — s> oo leads to the same result for the stationary phase condition and correlation of u(s) 
as the previous one in Eq. (p^. In section 4, we adopt this approach to examine the elastic 
response of a semiflexible chain under tension. 



2.5 Operator representations 

For periodic chains, we can exploit the analogy with quantum mechanics and use an op- 
erator representation of the free energy. This is especially useful as L ^ oo which allows the 
use of ground state dominance approximation. If we interpret p = Zp^u as the momentum 
operator and u as the position operator such that they satisfy the commutation relation 
[pj,Uk\ = —iSjk, the free energy of periodic chain can be written as 

J^[\] = -In tr e-^^[^l - LA (44) 

where the Hamiltonian 7i is defined by 

m = 1^ + Au^ (45) 

This problem is equivalent to a 3-dimensional quantum harmonic oscillator moving in imag- 
inary time s = —it; if we denote the energy eigenvalues of the Hamiltonian in Eq. (^) by 
En with n = (ni, 77-2, ^3) and rii = 0, 1, 2..., then 

00 

Ln[x] ^ e"-^'^" 



tre-^^W = J2 



ni,n2,n3=0 

^g-LO{n+i) 


(2^h^h(pzy; ■ 

This calculation gives the same stationary phase condition as before. 



(46) 



15 



This approach is especially useful as L ^ oo. For large L, we can use the following 
formula 

^-LH ^ ^-LEo J|o^^o| + c(e-^(^^-^o))] (47) 

where Ei = |f2 and Ei = The ground state denoted by |0) corresponds to the lowest 
eigenvalue Eq of H and is assumed to be unique and separated by a gap from Ei. As 
L — > oo, this expression is dominated by the ground state. The term tr exp(— LTi) can be 
easily computed to yield the stationary phase condition, i.e., = |- Introducing U(s), 
the interaction representation of the operator u, 

U(s) = e"^^ u e^^ (48) 

we can express the correlation of u as 

, , , , tr[U(s) ■ U(s')ce-^^] , ^ 

{u{s) ■ u{s')) = y ^, 0<s<s'<L. (49) 

tr e~^'^ 

As L ^ oo, the ground state is dominant in the above expression contribution: 

(u(.)-u(.')) = (o|ru(.) -11(5010) (50) 

where T is a time ordering operator. In general, the correlation functions of the statistical 
analog correspond with the time-ordered products of the corresponding quantum fields. If 
we use the basis |n) in which 7i is diagonal, the correlation can be further simplified to yield 



(u(s') ■ u{s)) = (0|ue-l^'-^l^ u el^'-^l^|0) 
= ^|(0|u|n)|V(^"-^o)l^'-^l 



^-{Ei-Eo)\s'-s\ 



= e-l^'-^l/^o. (51) 

As we have seen, the operator representation of the free energy are especially useful for a 
very long chain. In general, we can make ground state dominance approximation for a long 
chain as long as the Hamiltonian of the statistical analog Ti. is hermitian bounded below and 
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has a discrete energy spectrum. In this case, the free energy per unit length is approximately 



to calculation of transition amplitudes (0|u|n). For an ideal semiflexible chain, u connects 
(0| with the first excited state with no n- vanishing transition amplitude, resulting in the cor- 



3 END-TO-END DISTRIBUTION FUNCTION 

There have been several studies over the years which have focused on the distribution of 
end-to-end distance of semifiexible chains. Most of these studies have started by considering 
the chains near the rod limit and have computed corrections in powers of where t is 
the ratio of the bare persistence length to the contour length. In general the calculations 
have been done only to low orders in t^^. These calculations are quite involved, and more 
importantly, they do not provide reliable results in the interesting cases, i.e., when t is on 
the order of unity. In a very recent paper Wilhelm and Frey have reported analytic and 
numerical (Monte Carlo) calculations for the radial distribution function for a range of values 
t. The analytical expressions are given in terms of an infinite series involving the second 
Hermite polynomials for chains in three dimensions. The meanfield like theory presented 
in the previous section can be used to derive a very simple expression for the distribution 
function of end-to-end distance. 

The distribution function of the end-to-end distance R is 



equal to the ground state energy of H. The calculation of the correlation function reduces 





(52) 



where the average is evaluated as 




/D[u(s)]...Vi>MF[u(g)] 

JV[uis)] v|/mf[u(s)] 



(53) 
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The weight ^[u(s)] is (see Eq. (|)) 



\E'[u(s)] oc exp 



(54) 



0<s<L 



Following the previous section, at the level of the stationary phase approximation, the weight 
can be replaced by 



^mf[u(s)] cx exp '^^i~^y~^Jo dsu2(^) --^(uo + ui) 



(55) 



where uq = u(0) and u^, = u(L). The parameters A and S, which are used to enforce 
the constraint u^(s) = 1 by a global constraint (u^(s)) = 1, are determined variationally. 
From the result in sec. 2.2, it is clear that due to chain fluctuations at the ends the global 
constraint (u^(s)) cannot be imposed using only one parameter as suggested by others p7 |. 
The variational solutions for A and 6 depend on the problem under consideration. For the 
calculation of G{R] L) the optimal values of A and 6 are different from that obtained in 
section 2. 

We calculate G{R;L) by replacing the true weight in Eq. (15^) by \E'mf[u(s)] given in 
Eq. (|55D. The equation for G{R; L) with the weight given by \1/mf[u(s)] is 

rioo 



G{R-L)=U-' f^^j-yj dAd5 j V[u{s)] e-'^-(i^-/o'd-(^)) vl>j^p[u(.)] 



(56) 



where is an appropriate normalization constant. The functional integral over u(s) in 
Eq. (^) is done by replacing u(s) by u — ^. The resulting path integral corresponds to a 
harmonic oscillator that makes a transition from Uq — ^ to u/^ — ^ in "time" L. Using the 
standard result for harmonic oscillator propagator the distribution function G[R] L) becomes 



G{R]L) = JdxJdS exp{-^[A,(5]} 



(57) 



where 



In 



LA + 2(5^ 



\ , /sinhf]L\ , , ^, 
)+ln(^^)+ln(a/3) 



4A2 



(58) 



with 



a = + Lp\l cotnl 



fnL\ 
I 1 



and 



2A 



(59) 
(60) 

(61) 



We evaluate the integrals over A and 5 (see Eq. (pT])) by the stationary phase approximation. 
This approximation replaces the constraint u^(s) = 1 by a global constraint (u^(s)) = 
1. If the global constraint is enforced using only one variational parameter, as suggested 
elsewhere l^^, then one simply obtain the incorrect Gaussian expression for G{R;L). The 
presence of the parameter 6 accounts for the suppression of the fluctuations of the ends of 
the chain. The stationarity condition 

d 



gives, after some algebra. 



I XL 



J^[A, 6] = 



1 



1 _ ^ 



(62) 



(63) 



and we find that 6 does not show any significant variation with R. Using the stationary values 
of A and 6 the distribution of the end-to-end distance for the semifiexible chain becomes 



G{R; L) = const - 



1 



1 - ^r/' 



exp 



9L 



1 



% (1 - f 



(64) 



We showed in the previous section that within the stationary phase approach, which enforces 
thae constraint only globally, the persistence length is reduced to /q = instead of Ip. If we 
let t = ^ the radial probability density in three dimensions for semifiexible chains is given 



by the simple expression, namely. 



P{r-t) = AnC- 



i)9/2 



exp 



3t 1 
T(l -r2 



(65) 
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where r = ^. The normahzation constant C which is determined using the condition 

drP{r;t) = 1 (66) 

Jo 



and is given by 

1 

7r3/2e-"a-3/2(lT3a-i + fa' 
where a = St/A. 

The distribution function, P(r; t), goes as as r — and vanishes at r = 1. The peak 
of the distribution function occurs at 



rmax = y ^ (68) 

where ?7 = | — |t. In Fig. (1) we plot P(r; t) for the five values of t for which Wilhelm and 



Frey |^ have presented simulation data. We find that our simple expression in Eq. 
almost quantitatively reproduces the data with the maximum deviation of about 10% at the 
peak for t = 0.5. 

We have also calculated the first two moments of the distribution function and find that 
they reproduce the exact results in both the random coil limit and the rod limit. The distri- 
bution function P(r; t) also has the correct limiting behavior (~ 6{1 — r)) as t ^ 0. Thus the 
meanfield variational approach yields a simple expression for P(r; t) that is in quantitative 
agreement with simulation and hence can form the basis to analyze experiments. 



4 SEMIFLEXIBLE CHAINS UNDER TENSION 

The meanfield theory described in the previous section is especially useful when one 
encounters semifiexible chains in the presence of fields for which tractable calculations are 
difficult. In this section we use the variational theory to investigate the external field on 
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the conformational properties of a semiflexible chain. Even though the analogous problems 
in flexible chains was extensively studied, the corresponding problems in the semiflexible 
chains has only recently attracted much attention [27-32]. The earliest theoretical paper 
dealing with this problem was initiated by Fixman and Kovacs [^. These authors used 
a modified version of the Gaussian model for stiff chains and provided expressions for the 
stretching as a function of applied force. Their treatment is only valid when the applied 
force is small and significant deviation from these predictions are observed at sufficiently 
large values of the external force. Marko and Siggia [31-32] have calculated the extension 
as a function of force for wormlike chains and found that their results fit the experimental 
data very well. Some aspects of this theory have also been considered by Odijk [^ who 
also discusses the competition between entropically dominated effects and elasticity effects. 
The meanfield type approach adopted here is, we believe, more general than those adopted 
in the literature. 

Gaussian chains can be arbitrarily extended under tension / as implied by Hook's law 
of elastic response, i.e., R ~ Llf; real chains cannot be extended beyond the contour length 
L. As the magnitude of tension exerted on both ends of a semiflexible chain increases one 
expects an interesting crossover to occur from Hook's limit to fully extended limit {R = L). 
Thus the problem of stretched semiflexible chains entails a competition between entropy 
dominated effects and ordering due to the external field f . We will show that the meanfield 
treatment quantitatively reproduces experiments on DNA [^ thus providing further evidence 
for the validity of meanfield approach. 



4.1 Stationary phase condition 

For mathematical convenience, we will adopt translationally symmetric chain model in- 
troduced in Eq. (pS]). Each segment in a semiflexible chain is assumed to be stretched by 
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tangential force f(s) which tends to suppress chain fluctuations. The weight for the case of 
semiflexible chains in the presence of the stretching is 

^[u(s), \{s)] oc exp [-i r r dsds'u(s)Q(s, s')u(s') + u(s) ■ f (s)dsl (69) 
I Jo Jo Jo J 

where Q is the same operator as introduced in Eq. (|39|) , i.e., 

6{s'-s). (70) 



r 

Q{s,s')= -/p(-)' + 2A 



The last term in the exponent is the energy penalty for chain conformations which are not 
parallel to f(s). The free energy can be written as 

exp(-F/A;BT) oc H V[X{s)] [ V[u{s)] e^o ^u{s) , X{s)] 

J —too J 

(X £2^[A(s)]exp{-^[A(s),f(s)]} (71) 

where JF is the free energy functional given by 

^[A(s), f (s)] = -In J P[u(s)]^[u(s), A(s)] - A(s)ds + const.. (72) 

Integration with respect to u(s) yields 

^[A(s),h(s)] = ItrlnQ -l££i{s)Q-\s,sy{s') 

- X(s)ds + const.. (73) 
Jo 

Note that jF[A(s), f (s)] is the generating functional for the connected correlation function 
defined by 

6 8T 



(u(.) ■ u(.')), 



(5f(s) 5f(s') 
d 

■ (u(s')) 



5f(.) 

= (u(.).u(s'))-(u(s)).(u(/)). (74) 

As f goes to zero, (u(s)) vanishes. In this case, the distinction between the connected 
correlation function and the correlation function disappears. Performing the functional dif- 
ferentiations in Eq. (46), we get 

(u(.)-u(.'))e = 3g-^(s,s') (75) 
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and 

M.)) = -^ = I Q-\s,s'W). (76) 

Integrations of each term in Eq. (^) with respect to s and s' lead to the following expression 
for A/?2 

Ai?2 = {R^)-{Rf 

= ["^ r{u{s) ■ u{s'))dsds' - r{u{s)) ■ {u{s'))dsds' 
Jo Jo Jo Jo 

= JoJoWs)' Wf'"^''- ^^'^ 
All relevant quantities can be thus expressed in terms of the free energy functional jF[A(s), f (s)]. 

Following the general formalism described in section 2, the case of ideal semiflexible 
chains, the free energy is approximated by the minimum value of the corresponding free 
energy functional which occurs along the stationary phase trajectory A. The stationary 
phase condition is obtained by requiring j^^jJ^ = 0. To this end, we first note that 

^ -trlng = ^E(lnQ)(^''^') 



and 



5X{s) ^ 5\{s) ^ 

= 2Q~\s,s). (7J 



^ Q-\s',s") = 2Q-\s',s)Q-\s,s"). (79) 



5X{s 

The stationary phase condition leads to 

1 = £ £ ds'ds"iis')Q-\s',s)Q~\s,s")fis"). (80) 

The stationary phase value A(s) thus depends on the value of the external force f(s). The 
simplest but nontrivial case corresponds to the case of f(s) = f = const. It is the case of 
stretching of semiflexible chains under constant value of f that is appropriate to the recent 
experiments on DNA [^]. If f is uniform, then so is A. In this case the stationary phase 
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condition reduces to 



2\ds'^ "T ^/ s,s' 



IE 



Ip (27rn)2 



oo AL + L 



^ 4A2 
4A2' 



^1) 



With the aid of Eq. (|33|) , this can be written in a closed form 



i = jgcoth(inL) + ^ 



i2) 



where = As L — > oo, this equation can be further simphfied to yield 



1 - 



^3) 



As / — + 0, this equation reduces to the earlier one in Eq. ( pTj ) as expected. 



4.2 Conformation of a semifiexible chain under tension 



In the case of uniform f , we obtain a simple expression for the correlation function: 



(u(.)-u(.')) = 3Q-\s,s') + ^^ 



3 g exp[z^|s' - s\ 



2 X r , ip (27rn)2 + 



-oo AL + 



2 L 



4A2 



/pA 



2 cosh[(L-2|s'-s|)fi/2] 

sinh(fiL/2) ^ 4A2' 

where Vt = \J^- To derive the last step from the previous one, we used Eq. 
the above equation can be further simplified to yield 



(84) 

As L — > oo, 



{n[s) ■ u{s')) 



"V /nA 



+ 4A2- 



^5) 



If we set |s' — s| to zero, this equation reduces to the stationary phase condition. Thus 
the free energy minimization condition in this case also amounts to requiring u^is) = 1. As 
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f — i> 0, the coefficient in the exponential function becomes 1, resulting in the same expression 
as in Eq. (|5TD . 

Using the correlation function in Eq. (^), the mean squared internal distance of the 
chain under tension can be obtained 



(|r(.')-r(.)p) = r'r'(u(5i)-u(.2))d.id.2 

J s J s 



'A 



l\s' -s\~—(l-e-^^''-'^' 



+ ^^is'-sr. (86) 



4A 



We can easily see that, for small \s' — s\, the ffist two term on the R.H.S. of Eq. (|86|) is 
dominant. As |s' — s| becomes larger, the last term is more important. At large length 
scales, the chain conformation is thus mainly determined by the orienting field. We can now 
introduce a crossover length Sf at which these two length scales coincide. If the crossover 
occurs at Sf which is somewhat larger than then we can get a simple expression for Sf] 
this can be obtained by balancing the ffist term with the last one in R.H.S of Eq. (186|) 



2 / 4 



Sf = hll^(i-.^n]- (87) 



The condition, Sf ^ Q^^, is equivalent to ^ <^ \\fv^ ~ ^(1) f^v ^ 1- Under this 
condition, the stationary phase condition coincides with that for an ideal semiflexible chain. 
This results in 

Sf ~ (88) 



The longitudinal size corresponding to Sf is called a tensile screening length, ||3^ and is 



denoted hy S,f- By noting that parts of the chain within the length scale Sf resemble an 
ideal chain, we can have 

if ~ f-\ (89) 

The condition <^ 1 for the above expression to be valid in turn amounts to requiring 
//p -C 1 or equivalently ^/ ^ /p. It is interesting to note that results for Sf and ^/ given 
above are the same as those for a Gaussian chain under tension. This can be understood 
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as follows; the condition, ^/ ^ /p, for the above results to be valid, is equivalent to saying 
that the tensile screening length contains large number of chain segments of length /p. 
The chain stiffness becomes marginal in determining the shape of parts of chains within the 
length scale ^f. Since the effect of f is important only beyond this length scale, the chain 
stiffness is not "coupled" to the external field f . 

As ^ becomes much larger than f or equivalently flp oo and thus Sj approaches 
0. This implies that, at any length scale beyond Sj ~ 0, the chain conformation is governed 
by the interaction term — / f ■ u(s). As we will see below, as Ipf —>■ oo, the entropical contri- 
bution is not balanced by the energy penalty for crumpled conformation. That is the chain 
segments are aligned with f as expected. These quantitative considerations, which were 
already implicit in the calculations of Fixman and Kovacs p8|], can be made quantitative 
using our general formalism. 



4.3 Elastic response of a semifiexible chain 



Let us now consider the average elongation z defined by 



L f 
dsu(s) ■ - 
/ 



(90) 



This quantity, which is experimentally measureable, shows how the chain responses to the 
tangential field. The statistical average in Eq. ( pOD can be conveniently expressed in terms 
of the free energy functional JF (see Eq. (0)). For the case of uniform f, we obtain a simple 
expression for this: 



f 5i ^ Jo 

= f 



/ 5 \ 2 



ds 



2A(/) 



(91) 



where the argument in A(/) is introduced to emphasize the dependence of A on /. This 
equation along with the stationary phase condition Eq. (R3) provide the average elongation 
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with respect to /. To obtain analytic expressions, we first take two limiting cases of /p — > 
and / ^ oo. As /p — > or more precisely Ipf — > 0, the problem reduces to that of a 
stretched Gaussian chain as mentioned in sec. 4.1. In this case, we can get the following 
elastic response relation 

z = /f (92) 

where = 21qL = 2(|/p)L is the size of the corresponding ideal semiflexible chain. This is 
the elastic response relation of a Gaussian chain under tension [^. For the case of flp <C 1, 
the stationary phase condition can be expanded in terms of p. Up to p, we get 

m ^ A + ^ (93) 

where A is the stationary phase condition for / = 0. Thus the average elongation is given by 

z^f§[l-im (94) 

Note here that the control parameter in this expansion is flp (or flo). 

As /p increases the term in the stationary phase condition becomes important. If 
the value of /p reaches the same order of magnitude as that of we expect the effect of 
the orienting field to be important over all length scales. In other words, chain fluctuations 
become frozen as Ipf becomes large. As / — > oo, the stationary phase condition has a 
solution at A^^ = 0, resulting in A = |/. The average elongation is thus given hy z = L. 
Chain fluctuation are totally suppressed in this case. For finite but much larger / than /p, 
A(/) can be approximated as 

M/) = i/(l + !^)- (95) 
Accordingly z for finite / is smaller than L, which is due to chain fluctuations. 



1 



z^L[l-i^-). (96) 

Similar result with slightly larger numerical factor was obtained by T. Odijk considering 
a semiflexible chain near the rod limit. 
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Since the chain under strong elongation as described by Eq. ( p6| ) is nearly parallel to the 
orienting force f , we can approximate z as 

/ COS 6'/(s)ds 

JO 

L 



\-WAs))\ (97) 

where ^(s) is the angle between f and u(s). For a uniform field, we expect (dfis)) not to 
depend on s for a long chain, i.e., {Oji^s)) = (Oj). Comparing Eq. (|96| ) and (p7|) , we can 
obtain the expression for this quantity: 



m - i^jj. (98) 

This result again differs by numerical factor from the corresponding expression for (Oj) 
derived by Odijk . 



In the above analysis, we have seen that the stationary phase approach reproduces results 
consistent with literature in limiting cases. To test our theory further, we solved the station- 
ary phase condition numerically. To compare with experimental data obtained by Smith et 
al. on the stretching of DNA (see also Fig. 2 in Ref. |^2|]) Iq and L have been chosen to 
be 53.4nm and 32.80/im respectively. In Fig. 2, z/L is plotted against In /. Our numerical 
solution (continuous curve) is well in quantitative agreement with the experimental data 

In the previous analysis, we have seen that a strong orienting field tends to suppress 
chain fluctuations as is implied by the correlation and the elastic response relation. This is 
also manifested in the ratio Using Eq. ([77|), we get 

2 , 2 



4' 



—L - —d - e^^^l 



and thus 



(99) 



(100) 



28 



For a short chain, = 0{1) while it approaches zero as L — oo. It can be easily seen 

that value of this quantity crosses over at the contour length L Sf from fluctuation dom- 
inated limit {-^0^ = to fluctuation suppressed limit {^0^ ~ 0). As l^f oo, Sj 
approaches zero. Thus in this case, chain fluctuation becomes totally suppressed for a chain 
with contour length L ^ ~ 0. 



4.4 Operator representation 

Just as in the case of an ideal semiflexible chain, it is instructive to use operator repre- 
sentation of the free energy functional J^[A, f]. The corresponding Hamiltonian in this case 
describes a quantum harmonic oscillator in the linear force field described by f: 

7i:[A,f] = ?i:[A] + A7i: 

Since [p, f] = and thus the commutator [pj, Uk\ = —iSjk is invariant under u ^ u — the 
external field just shifts eigenvalues of 7i[A, f] by the same amount This results in 

This leads to the same stationary phase condition as in Eq. (0). 

Correlation function of u(s) can be also computed following the similar steps to that 
leading to Eq. (^Tj) 



(u(.)-u(.')) = E|(0|(u+^)|n) 



2 ^-{En-Eo)\s'-s\ 



P 

' -h - 

z 

3./^e-l^'-^l^ + -^. (103) 



|(0|u|n)|2e-(^-^»)l^'-^l + ^^2 

^ -W-s\n , _f 
^V^pA 4A2 
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The operator representation again produces the same result as the functional integral formal- 
ism. One of the advantages of using the operator representation is that we can use ground 
state dominance approximation provided there is a energy gap. 



5. SEMIFLEXIBLE CHAINS IN A NEMATIC FIELD 

As a final example of the utility of our approach we consider an semiflexible chain in 
a nematic environment. These calculations expose conditions under which the stationary 
phase approach is not successful. The interaction energy of a chain in a nematic field along 
n is assumed to be —g Jq[u{s) ■ n]^, where the coupling constant g measures the strength of 
nematic field. Thus conformations of nematic polymers either parallel or anti-parallel to the 
nematic field are equally energetically favorable. This is in contrast to the case of a stretched 



chain. The Hamiltonian of a single chain in the present case is given by 0, q6 



where u(s) is defined only on a unit sphere of |u(s)| = 1. At the meanfield level this takes 
the form 

= i / d^d^' "(^) '^(^') '9 f^^ds ulis) (105) 

/cb-Z jo jo jo 

such that the weight is given by 

v^MF[u(.),<7]cxexp|-^^} (106) 

where the direction of the nematic field is chosen to be parallel to z-axis. With the convention 
u_L = {ux,Uy,0), this can be rewritten as 



r rdsds'[ux(s) Q{s,s') ux(s') +M.(s) Qi(s,s') u,{s')\ (107) 

JO JO J 
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where 



and 



S{s' - s) 



-/p(|)V2(A-,) 



5(s'-s). 



(108) 



(109) 



The corresponding generating functional T[X, g,i{s)], up to an additive constant, now has 
the following form 



J^[X, = tr In Q - i dsds'f^{s)Q-\s, s')f±(s') 

Jo JO 

+ itr In - i dsds'Us)Qi\s, s')/.(s') - £ X{s)ds (110) 



from which we can obtain 



1 ^ 



- + (u^(.))-(u^(.')) 



(Ill) 



and 



I ^ ^-\s'-s\U 

V ^P(A - 9) 



- + («,(s))-(«,(s')) 



(112) 



with Q = and Jli = the last steps of above equations, we set f(s) to zero. 

The stationary phase value of A now satisfies the following equality 



which ensures 



(113) 



(114) 



It can be easily seen that the stationary phase value of A in this case is larger than that 
for the case of free chains. But (A — g) is smaller than that for the latter case. This implies 
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that the nematic field induces ordering along the field but decrease the persistence length 
perpendicular to it. That is 

(Rl) = n\n^{s) ■ u^{s'))dsds' 
Jo Jo 

< l{2loL) (115) 

and 

(Rl) = N\u,{s)-u,{s'))dsds' 
Jo Jo 

< W^loL) (116) 

where Iq — |/p is the persistence length of a free chain. The size of chains in nematic field 
does not grow as L as L ^ oo, as is the case for the stretched chain within the meanfield 
variational theory. To be more precise, let us consider the stationary phase solutions for 
two extreme cases. In the following derivations any quantity with (without) argument g 
corresponds to a nematic polymer (non-interacting polymer), 
(a) The weak nematic limit, g ^ 0. 

For small g, we can write X{g) ~ A(l + e). By assuming both e and g small, we can have 

A(^) ^ A(l + yio) (117) 

and thus 



n{g) = p\{g)/l^ ~ (1 + iglo) (118) 



and 



These lead to 



n'{g) = J2{\{g) - g)/l, ^ pX/l^ (1 - \gk). (119) 



(120) 
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and 

lffi)-^ + 59'o. (121) 

Thus the persistence length along the nematic field is increased by the factor of (1 + ^glo) 
which can be compared with (1 + ^glo) (in our notation) obtained by Warner at al. [^]. 
(b) The strong nematic limit, g oo. 

If we simply let g —>■ oo in the stationarity condition in Eq. (|113| ), then X{g) and thus 



fl{g) should diverge so that {X{g) — g) remain positive. This results in 

(ux(s)-ux(s')) = (122) 

and 

Ms) ■ u,{s')) = exp(-|s' - s|/2/p). (123) 

The diverging nematic field totally suppresses chain fiuctuations perpendicular to it but 
increases the persistence length parallel to it by the factor of 2. The increase of the persistence 
length along the nematic field by only numerical factor, however, is not correct. Several 
authors predicted an effective persistence which grows exponentially with J glp H) 12, 36|, 37 



Following the reference [^, this can be explained as follows; since the interaction energy 
~gu1 which is quadratic in Uz has two minima for = ±1, the "tunneling" probability from 
one minimum to the next is very small for large g. This results in the length scale which 
varies exponentially as \J'glp over which the chain is parallel to the nematic field. For more 
details see the reference [ p6| . 

The stationary phase free energy in Eq. ( |105| ) also contains the quadratic interaction 
term —gu^. This term if added to Xul becomes (A — g)u1 with positive (A — g) which has 
only one minimum at = 0. Thus the stationary phase approach does not account for the 
instanton solution, resulting in much smaller persistence length than the exact results. This 
problem illustrates the possible limitations of the stationary phase approach; if there is a 
possibility of broken symmetry in the problem then the meanfield variational theory should 



33 



be used with caution. 



6. CONCLUSIONS 



In this chapter we have described a simple meanfield variational approach to study a 
number of properties of intrinsically stiff chains which are appropriate models for a large 
class of biopolymers. The exact statistical mechanics of such systems are complicated by the 
fact that the constraint of the tangent vector being unity has to be enforced at all points 
along the contour. In the meanfield approach this local constraint, which is difficult to 
impose, is replaced by a global one. The global constraint ensures that the condition of the 
tangent vector being unity is obeyed on an average. We have described the calculation of 
the distribution of end-to-end distance and the elastic response of stiff chains under tension 
using this basic methodology. In the former example we find that the simple expression 
almost quantitatively fits the results of simulation. For the case of the stiff chain under 
tension we recover analytically all the known limits. Furthermore we obtain quantitative 
agreement with recent experiments on the stretching of DNA using the total contour length 
and the persistence length as adjustable parameters. The limitations of our approach become 
obvious in situations that involves broken symmetry such as the case of a stiff chain in a 
strong nematic potential. Nevertheless it is clear that the meanfield variational approach lays 
the foundation for systematic studies of more complicated problems involving semiflexible 
chains. 
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